Interaction induced fractional Bloch and tunneling oscillations 
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We study the dynamics of few interacting bosons in a one-dimensional lattice with dc bias. In the 
absence of interactions the system displays single particle Bloch oscillations. For strong interaction 
the Bloch oscillation regime reemerges with fractional Bloch periods which are inversely proportional 
to the number of bosons clustered into a bound state. The interaction strength is affecting the 
oscillation amplitude. Excellent agreement is found between numerical data and a composite particle 
dynamics approach. For specific values of the interaction strength a particle will tunnel from the 
interacting cloud to a well defined distant lattice location. 
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Bloch oscillations [l*] in dc biased lattices are due to 
wave interference and have been observed in a number 
of quite different physical systems: atomic oscillations 
in Bose-Einstein Condensates (BEC) [2|, light intensity 
oscillations in waveguide arrays J3|, and acoustic waves 
in layered and elastic structures y , among others. 

Quantum many body interactions can alter the above 
outcome. A mean field treatment will make the wave 
equations nonlinear and typically nonintegrable. For in- 
stance, for many atoms in a Bose-Einstein condensate, a 
mean field treatment leads to the Gross-Pitaevsky equa- 
tion for nonlinear waves. The main effect of nonlinear ity 
is to deteriorate Bloch oscillations, as recently studied 
experimentally [s] and theoretically [6|-|8|- 

In contrast, we will explore the fate of Bloch oscilla- 
tions for quantum interacting few-body systems. This is 
motivated by recent experimental advance [9] in moni- 
toring and manipulating few bosons in optical lattices. 
Few body quantum systems are expected to have finite 
eigenvalue spacings, consequent quasiperiodic temporal 
evolution and phase coherence. In a recent report on in- 
teracting electron dynamics spectral evidence for a Bloch 
frequency doubling was reported [lo|. On the other hand, 
it has been also recently argued that Bloch oscillations 
will be effectively destroyed for few interacting bosons 

ini- 

In the present paper we show that for strongly interact- 
ing bosons a coherent Bloch oscillation regime reemerges. 
If the bosons are clustered into an interacting cloud at 
time t = 0, the period of Bloch oscillations will be a frac- 
tion of the period of the noninteracting case, scaling as 
the inverse number of interacting particles (Fig(T]). The 
amplitude (spatial extent) of these fractional Bloch oscil- 
lations will decrease with increasing interaction strength. 
For specific values of the interaction, one of the particles 
will leave the interacting cloud and tunnel to a possi- 
bly distant and well defined site of the lattice. For few 
particles the dynamics is always quasiperiodic, and a de- 
coherence similar to the case of a mean field nonlinear 
equation [7] will not take place. 
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FIG. 1: Time evolution of the probability density function 
(PDF) Pj{t) for the interaction constant [7 = 3 and dc field 
E — 0.05 and different particle numbers initially occupying a 
single site at t = 0. (a) shows one particle Bloch oscillations 
with the conventional Bloch period 2ii/E^ while (b), (c) and 
(d) display two, three and four particle oscillations with the 
periods 2n/(2E), 27r/(3£;) and 2n/{AE), respectively. 



We consider the Bose-Hubbard model with a dc field: 
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where bf and bj are standard boson creation and anni- 
hilation operators at lattice site j; the hopping ti = 1; 
U and E are the interaction and dc field strengths, re- 
spectively. To study the dynamics of n particles we use 
the ortho normal basis of states |k) = /c2, ^n) = 
^ki^k2"'^kn 1^) where |0) is the zero particle vacuum state, 
and ki < k2... < kn are lattice site indices (for instance, 
in the case of two particles the state representation is 
mapped to the triangle). The eigenvectors \u) of Hamil- 
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tonian ([T]) with eigenvalues are then given by: 

W) = J2A'^\k), nW) = xM (2) 

k 

where the eigenvectors = (k|i^) and the time evolution 
of a wave function is given by 



{H*(0)). (3) 



We monitor the probability density function (PDF) 
Pj{t) = (^(t)|6^6j|^(t))/n, which can be also computed 
using the eigenvectors and eigenvalues: 
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In Figd] we show the evolution of Pj{t) for [7 = 3, 
E = 0.05 and n = 1,2,3,4 with initial state ki = k2 = 
... = kn = i.e. when all particles are launched on 
the same lattice site p. For n = 1 we observe the usual 
Bloch oscillations with period T = 27r/E (FigUJa) and 
below). Due to the small value of the amplitudes of 
oscillations are large. However, with increasing number 
of particles, we find that the oscillation period is reduced 
according to 27r/(n£^), and at the same time the ampli- 
tude of oscillations is also reduced. 

One particle case: For n = 1 the interaction term 
in ([T]) does not contribute. The eigenvalues Aj^ = Eu 
(with ly being an integer) form an equidistant spectrum 
which extends over the whole real axis - the Wannier- 
Stark ladder. The corresponding eigenfunctions obey the 
generalized translational invariance A^'^^ = A)^ f]\ and 
are given by the Bessel function Jk{x) of the first kind 
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Al = r,^J,_,{2/E). 



(5) 



All eigenvectors are spatially localized with an asymp- 
totic decay (l/E)^ /k\^ giving rise to the 
well-known localized Bloch oscillations with period Tb = 
27t/E. The localization volume £ of a single particle 
eigenstate characterizes its spatial extent. It follows 
C oc -[E ■ \nE]-^ for ^ and £ ^ 1 for E ^ oo 
0. For E = 0.05 the single particle oscillates with am- 
plitude of the order of 2C ^ 160 (FigUJa)). According to 
Eqs. (j4|) and (|5]) the probability density function is given 
by: 



(6) 



Two particle case (n=2): For [7 = the eigenfunc- 
tions of the Hamiltonian ([T]) are given by tensor products 
of the single particle eigenstates: 



ii<y. (7) 



FIG. 2: Eigenvalue spectrum for n = 2, ^ = 0.5 and different 
interaction constants U . The eigenvalues are displayed only 
for eigenvectors localized in the center of the lattice (we se- 
lect the 32 eigenstates which overlap most strongly with the 
center of the lattice) (a): [/ = 0, the spectrum is equidistant 
with spacing E and degenerate, (b): U — 2^ the degeneracy 
is lifted, (c): U — 15, the spectrum decomposes into two sub- 
spectra, with two different equidistant spacings - E and 2E. 
Graph (d) displays the eigenvalue spectrum of the 32 central 
eigenfunctions as a function of U . 



The corresponding eigenvalues form an equidistant spec- 
trum which is highly degenerate: 



H\pi,v) = {^^v)E\^,v) 



(8) 



For the above initial condition ki = k2 = p the expression 
for the PDF (|6]) is still valid (actually it is for any number 
of noninteracting particles), with the same period 27r/E 
of Bloch oscillations as in the single particle case. 

For nonvanishing interaction the degeneracy of the 
spectrum is lifted, and the eigenvalues of overlapping 
states are not equidistant any more (Fig 12]). Therefore we 
observe quasiperiodic oscillations which however are still 
localizing the particles. For even larger values of U the 
basis states with two particles on the same site will shift 
their energies by U exceeding the hopping 2ti. There- 
fore for U > 2ti the spectrum will be decomposed into 
two nonoverlapping parts - a noninteracting one which 
excludes double occupancy and has equidistant spacing 
E, and an interacting part which is characterized by al- 
most complete double occupancy and has corresponding 
equidistant spacing 2E^ which is the cost of moving two 
particles from a given site to a neighboring site. Some ini- 
tial state can overlap strongly with eigenstates from one 
or the other part of the spectrum, and therefore result 
in different Bloch periods. In particular, when launch- 
ing both particles on the same site, one strongly overlaps 
with the interacting part of the spectrum and observes a 
fractional Bloch period 27t/{2E). 
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FIG. 3: Upper plots: PDF for E = 0.1, n = 2, single site ini- 
tial occupancy and different interaction constants. For U — 
we find single particle Bloch oscillations. For U — A fractional 
Bloch oscillations take place, in agreement with Lower 
plots: probability density of the evolved wave function (darker 
regions correspond to larger probablilities) after one half of 
the respective Bloch period. For [7 = the two particles are 
with equal probability close to each other and at maximal 
separation. For [7 = 4 the two particles avoid separation and 
form a composite particle which coherently oscillates in the 
lattice. In lower graphs we use triangle k < m mapping for 
indistinguishable two particle state representation (index m 
increases from the right to the left). 



with fractional period 27r/(2£^) as observed in Figdjb). 
The amplitude of the oscillations is reduced with 
increasing U since the hopping constant t2 is reduced 
(FiglSl). For U = 3 it follows t2 = 0.5, and together with 
the doubled Bloch frequency the localization volume 
should be reduced by a factor of 4 as compared to the 
single particle case. This is precisely what we find when 
comparing Fig(Ha,b): for n=l the amplitude is 160 
sites, while for n = 2 it is 40 sites. In the lower plots 
in FiglH we show the probability density of the wave 
functions |(^(t)|k)p after one half of the respective 
Bloch period in the space of the two particle coordinates 
with ki = k and k2 = m. For U = both particles 
are with high probability at a large distance from each 
other. Therefore the density is large not only for k = m 
(the two particles are at the same site), but also for 
/c = 5, m = 85 (the two particles are at maximum 
distance). However, for U = 4 we find that the two 
particles, which initially occupy the site p = 45, do not 
separate, and the density is large only along the diagonal 
k = m with 35 < A: < 55. (For U = the localization 
volume is ~ 20.) Therefore, the two particles indeed 
form a composite state and travel together. 

n particle case: We proceed similar to the case n = 2 
and estimate perturbatively the effective hopping con- 
stant for a composite particle of n bosons. For that we 
use the calculated width of the n-particle bound state 
band for £^ = [14]. In leading order of 1/U it reads 
0: 



In order to calculate the amplitude of these fractional 
Bloch oscillations, we note that for £^ = there ex- 
ists a two-particle bound state band of extended states 
with band width ^/U'^ + 16 - U [14]. For large U the 
bound states are again almost completely described by 
double occupancy. Therefore we can construct an effec- 
tive Hamiltonian for a composite particle of two bound 
bosons: 



+ 1^3 



2EjR+R, 



(9) 
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For n = 2 the above expression gives t2 — 2//7, the first 
expansion term of the exact relation for two bosons (p^Oj) . 
The corresponding composite particle Hamiltonian 



The PDF is given by 



R 



nEjR\Rj . (13) 



where and Rj are creation and annihilation operators 
at lattice site j of the composite particle (two bosons on 
the same site) with the effective hopping 



t2 = 



VU^ + 16-/7 



The corresponding PDF is given by 



p,{t) = J2A;A^^A-<A>;. 



(10) 



(11) 



The composite particle eigenvectors = 
Ju-p{2t2/{2E)) are again expressed through Bessel 
functions, but with a modified argument as compared 
to the single particle case. Bloch oscillations will evolve 



(14) 



and the composite particle eigenvectors A^ = 
Jiy-p{2tn/{nE)). Bloch oscillations will evolve with frac- 
tional period 27T/{nE) as observed in Fig(Hc,d). The 
amplitude of the oscillations is reduced with increasing 
U since the hopping constant tn is reduced. For U = 3 
and n = 3 it follows ts = 0.17, and for n = 4 we have 
^4 = 0.01. This leads to a reduction factor 18 and 400 
respectively as compared to the single particle ampli- 
tude and yields amplitudes of the order of 9 and 0.5 
respectively, which is in good agreement with the numer- 
ically observed amplitudes (10 and 2 sites respectively) 
in FiglUcd). 
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FIG. 4: Time averaged and normalized localization volume L 
of the wavepacket which emerge from two initial distributions 
as a function of [/ for ^ = 5. Red (grey) curve: two particles 
are launched on the same site. Orange (light grey) curve: 
two particles are launched on adjacent sites. Inset: PDF for 
U = 19.79, with clearly observed tunneling oscillations. 



Tunneling oscillations: For n = 1 the amplitude of 
Bloch oscillations is less than one site if £^ > 10 [7]. Thus, 
for n > 2 and increasing values of [/, the amplitude of 
fractional Bloch oscillations will be less than one site if 
EU^~^{n — 1)! > 10. Then, n particles launched on the 
same lattice site p will be localized on that site for all 
times. The energy of that state will be n((n — l)/7/2 + 
pE). If however one particle will be moved to a different 
location with site then the energy would change to 
(n - l)((n - 2)U/2^pE) + qE. For specific values of U 
these two energies will be equal: 



(n-l)U = dE , d = q-p. 



(15) 



In such a case, one particle will leave the interacting cloud 
at site p and tunnel to site q at distance d from the cloud, 
then tunnel back and so on, following effective Rabi os- 
cillation scenario between the states p^p) and |p, q). This 
process will appear as an asymmetric oscillation of a frac- 
tion of the cloud either up or down the field gradient (de- 
pending on the sign of U). We calculate the tunneling 
splitting of these two states using higher order perturba- 
tion theory, for an example see Ref.[15]. The tunneling 
time is then obtained as 



E^-\d-l)\ . 



(16) 



In order to observe these tunneling oscillations, we 
compute the time averaged second moment m2 = 

EPPjit) - {EjjPj{t)y of the PDF P. Then an ef- 
fective time-averaged volume of the interacting cloud is 
taken to be L = ^/12rn2 + 1 . We launch n = 2 parti- 
cles at site p = 40 and plot the ratio L(U)/L(U = 0) 
in FigH] (blue solid line). We find pronounced peaks 
at = E^ 2E^ 3E^ AE which become sharper and higher 
with increasing value ofU. As a comparison we also com- 
pute the same ratio for the initial condition when both 



particles occupy neighbouring sites (dashed red line), for 
which the resonant structures are absent. According 
to the above, the resonant structures correspond to a 
tunneling of one of the particles to a site at distance 
d = 1,2,3,4. The width of the peaks is inversely pro- 
portional to the tunneling time Ttun^ and the height in- 
creases linearly with the tunneling distance d. In the 
inset in FiglH we plot the time evolution of the PDF Pj 
for U = 19.79. We observe a clear tunneling process from 
site p = 40 to site g = 44. The numerically observed tun- 
neling time is approximately 1730 time units, while our 
above prediction ([16]) yields rtun ^ 1666, in very good 
agreement with the observations. 

Conclusions. The above findings can be useful for 
control of the dynamics of interacting particles. They can 
be also used as a testbed of whether experimental studies 
deal with quantum many body states. One such testbed 
is the observation of fractional Bloch oscillations, another 
one is the resonant tunneling of a particle from an inter- 
acting cloud. An intriguing question is the way these 
quantum coherent phenomena will disappear in the limit 
of many particles, where classical nonlinear and noninte- 
grable wave mechanics are expected to take over. 
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